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ABSTRACT 

The Decadal IRAC Bootes Survey (DIBS) is a mid-IR variability survey of the ^9 sq. deg. of the 
NDWFS Bootes Field and extends the time baseline of its predecessor, the Spitzer Deep, Wide-Field 
Survey (SDWFS), from 4 to 10 years. The Spitzer Space Telescope visited the field five times between 
2004 and 2014 at 3.6 and 4.5 microns. We provide the difference image analysis photometry for a half a 
million mostly extragalactic sources. In mid-IR color-color plane, sources with quasar colors constitute 
the largest variability class (75%), 16% of the variable objects have stellar colors and the remaining 9% 
have the colors of galaxies. Adding the fifth epoch doubles the number of variable AGNs for the same 
false positive rates as in SDWFS, or increases the number of sources by 20% while decreasing the false 
positive rates by factors of 2-3 for the same variability amplitude. We quantify the ensemble mid-IR 
variability of ^1500 spectroscopically confirmed AGNs using single power-law structure functions, 
which we find to be steeper (index 7 « 0.45) than in the optical (7 ~ 0.3), leading to much lower 
amplitudes at short time-lags. This provides evidence for large emission regions, smoothing out any 
fast UV/optical variations, as the origin of infrared quasar variability. The mid-IR AGN structure 
function slope 7 seems to be uncorrelated with both the luminosity and rest-frame wavelength, while 
the amplitude shows an anti-correlation with the luminosity and a correlation with the rest-frame 
wavelength. 

Subject headings: cosmology: observations — galaxies: active — quasars: general — infrared: galaxies 


1. INTRODUCTION 

There is increasing evidence that understanding the 
evolution of active galactic nucle0 (AGNs) is cru¬ 
cial for u nderstanding galaxy formation and evolu¬ 
tion (e.g.. iHonkins et al.l 120061 iDi Matteo et al.l 120081 
iKormendv fc Holl2013ll . This has driven a rapid growth 
in attempts to comprehensively characterize the proper¬ 
ties and demographics of AGNs. The Spitzer Space Tele¬ 
scope has played a significant role in this effort through 
the d evelopment of mid-IR AGN selec t ion methods 
(e.g., Lacv et id] 12004 iStern et al.l 120051 iHickox et al.l 


ies have expanded further based on the all-sky Wide- 




Donley et_ai.r 20081 iKozlowski fc KochanekI 120091 


Donley et al.l 1201^ and detailed studies of AGN spec¬ 
tral e nergy distributions (SEDs; vF,j) in the near/far-IR 
f e.g.. iRo wan-Robinson et Ml 120051 iRichards et al.l 12 0061 


Donley et al.l 1200^ Assef et al.l 2 0101 iKartaltene et alJ 

2 OIOI iPark et al.ll2010l iGhung et al.ll2014ll . These stud- 
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2010(1 mid-IR survey data (e.{ 

Ilchikawa et al. 

2012; 

Stern et al.l 120121 iTrichas et al.l 

120121: Assef et al.l 

2013; 

Yan et al. 20131 Mateos et al.ll2013HSatvaDal et al. 

2014; 

Assef et al.ll2015H. 


There is another dimension in which Spitzer can sig¬ 
nificantly expand its studies of AGN demographics and 
physics, namely the time domain. AGNs are vari¬ 
able sources at all wavel e ngths (e.g., iFalcke et al.|[2Mll: 
Vanden Berk et al.l 12004 iSoldi et 20081: ' Abdo et al 


2QQ9I: iMacLeod et al.l 120101 IKozlowski et al.l l201d.__ 
Denney et al.ll2010l: iGrier et al.ll2Q12[l . Their SEDs have 


s 


a minimum near a rest wavelength of I /rm and rise to¬ 
ward s both the “blue” and the “red” fe.s.. lElvis et al.l 
11994 lAssef et all 12010(1 . The blue, optical side corre¬ 
sponds to emission from the accretion disk surrounding 
a super massive black hole at the center of an AGN, while 
the red mid-IR side is due to emission by hot dust (e.g. , 
iShieldsI 119781 lAntonuc^ 119931: lUrrv fc Padovanil 11995(1 
as shown by near- lR reverberation mapping studies of 
nearby AGNs (e.g. iSugaiiuma et al.ll2004 Ikoshida et al ' 


20l4lShaDDee et al.ll20l4lVazauez et al.H2015HJun et~ 


2015(1 but also by direct observations of the torus 
(jAsmus et al.l 12014 1. Depending on redshift, the 3.6- 
4.5 pm emission probes from near the SED minimum 
(z > 3) into the dust-dominated regime (z < 1). 

We know a great deal about how quasars vary, but 
we know relatively little about the origins of that vari¬ 
ability. Possible origins include accretion disk instabili¬ 
ties, surface temperature fluc tuations, or variable heat¬ 
ing from coro i ial X- r avs fe.g.. iShakura fc SunvaevlIlQT^ 
Rokaki et al.l 119931 : iGhen fc Taaml 119951 : iRuan et al.l 
2014(1 . The introduction of stochastic process mod- 
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els, particularly the damped random walk (DRW; 


Kellv et al.l 120091: IMacLeod et al.l 120101: IKozlowski et al. 

2 OIOI: iButler & BloomI 20111: IMacLeod et al. 

2011 ; 

Zu et al. 

120111 IMacLeod et al.l 120121: IR.uan et al.l 

2012 ; 

Zu et al. 

120131). has provided a means to quantitativelv 


describe the variability of increasingly improving (in 
length and cadence) qu asar light curves from optical 
surveys (such as SP SS, iMacLeod et al.l 120101 OGLE . 
iKozlowski et al.ll201ClL La Silla-Q uest. lCartier et al.ll2015L 
or CRTS. iGraham et alJ l2015h or reverbe r ation map¬ 
ping campaigns (e.g., iDemnev et al.l 120101 : iGrier et ^ 
1201^ iDe Rosa et alJ [^15[ 1. The DRW model is char¬ 
acterized by two parameters, the time scale and am¬ 
plitude, where the time scale has been tied to eithe r 
disk orbital or thermal time scales (|Kellv et alJl200^ . 
While it well-describe s quasar light curves with a typ¬ 
ical cadence of days (jZu et al.l 1201,111 . there appears to 
be a break in the power spectrum distribution at short 
times scales observed for highly sampled Kep ler AGNs 
(jMushotzkv et alJ 120111: iKasliwal et al.l 1201^ . Sparser 
light curves for hundreds or thousands of AGNs are 
typically anal yzed with an ensemble structure function 
analysis ('e.g.. iHa wkinsI [^002t iVaii den Berk et al.l 12004 

Ide Vries et alJ[^05HSchmidt et al.ll2010ll . commonly de¬ 

scribed by a single power-law function with the slope in¬ 
dex 7 and amplitude Sq at a fixed time sc ale, although on 
long time scales there is clearly a break (IMacLeod et alJ 

[ml. 


In IKozlowski et ^ ()2010b[ l. we used a structure func¬ 
tion analysis to study the ensemble quasar variability 
of the largest repeatedly observed mid-IR AGN sample 
from the Spitzer Deep, Wide-Field Survey. The sur¬ 
vey observed the NDWFS Bootes Field ((RA, Deck) = 
(1430, -1-34:05)) four times between 2004 and 2008 
IT lAshbv et al.l[200l . It is a unique 9 deg^ area of the 
sky, far from the Galactic plane (^67 deg), that has 
been observed at many different wavelengths to provide 
a good census of extragalactic o bjects. The AGN and 
Galaxy Evolution Survey ((AGES: [Kochanek et al.ir2012tj 
provides redshifts for over 23,500 objects. The NO AO 
Deep Wide-Eield Survey fNDWFS: lJannuzi fc Deviri999ll 
provides deep optical images and photometry in the 
B^RI bands, while the Spitzer D eep, Wide-Field Sur¬ 
vey (SDWFS; lAshbv et al.l 12009( 1 provides magnitudes 
in four mid-IR bands at 3.6, 4.5, 5.8, and 8.0 mi¬ 
crons. It was a lso observed at X-r ay (XBootes; 


Murray et al.l l200l . at UV by GALEX (iMartin et alJ 
200511 at 24 nm (MAG ES survey), at far-I R by Her- 


schel Spa ce Observatory (iPilbratt et al.ll2010ll. and radio 
(FIR ST; iBecker et al.l 119951 and WSRT: Ide Vries et al.l 
l 2002 tl wavelengths. 

In this paper, we extend our previous analysis using the 
Decadal IRAC Bootes Survey (DIBS), which stretched 
the time baseline of the survey to 10 years by adding a 
fifth epoch. This expansion enables two improvements 
over our earlier study. First, the longer temporal base¬ 
line provided by the fifth epoch significantly improves 
the structure function slope measurement, our ability to 
characterize the differences between optical and mid-IR 
variability, and to subdivide the sample based on other 
criteria (e.g.. X-ray or radio flux). The ensemble quasar 
variability in mid-IR can be described by a single power- 
law structure function with the index 7 « 0.5, hence in¬ 


creasing the time baseline from 4 (SDWFS) to 10 years 
(DIBS), typically increases quasar variability by 50%. 
Second, existing mid-IR AGN samples are strongly bi¬ 
ased against low luminosity AGN and variability pro¬ 
vides a means of removing this bias to obtain a more 
complete census - the addition of the fifth epoch should 
improve the completeness and/or purity of the variability 
selection. Please note, however, that the typical rms of 
mid-IR quasar variability is ~0.1 mag, which is compara¬ 
ble to the uncertainties provided by DIBS or SDWFS for 
typical sources, hence we detect only significantly vari¬ 
able or brighter quasars. 

The paper is organized as follows. In Section [51 we 
present the available data and data analysis. In Sec¬ 
tion |3l we introduce the variability measures, while in 
Section [H we discuss the results of the quasar selection 
methods. The structure function analysis is presented in 
Section m and the paper is summarized in Section |6l 

2. DATA ANALYSIS 

A detailed desc ription of SDWFS is presented in 
lAshbv et al.l (120091 1. The first epoch was taken on 2004 
January 10-14 as the IRAC Shallow Survey (?). Then, 
we added three more epochs on 2007 August 8-13, 2008 
February 2-6, and 2 008 March 6-10 as the Spitzer Deep, 
Wide-Field Survey (|Ashbv et al.ll2009t D. Stern PI, PID 
40839). The new, fifth epoch was taken on 2014 April 23- 
29 as the Decadal IRAC Bootes Survey (DIBS; M.L.N. 
Ashby PI, PID 10088). The intervals between any two 
epochs span the range from 1 month to 10 years. 

Each epoc h consists of 7 ^ 20, 000 individual 256 x 256 
pixel IRAC (iFazio et al.l 1200411 frames that were com¬ 
bined into a single mosaic. We have re-reduced the 
original SDWFS images and created new mosaics, each 
approximately 19000 x 23000 pixels, with a pixel scale 
of 0.60 arcsec. These mosaics cover the same area of 
the sky, but at a different pixel resolution from SDWFS 
(13600 X16700 pixels and 0.84 arcse c). We used the sain e 
photometric catalogs as provided in lAshbv et al.l (1200911 . 
Due to the slightly different mosaicing of the fifth epoch 
and a slightly modified image subtraction procedure (de¬ 
scribed below), we report the photometric measurements 
for a total of 514975 objects. 

The photometry of our sour ces was done using dif - 
ference imaging methods (e.g., lAlard fc Luotopl 19^, 
where we use the implementation of iWozniakI ( 200011 . 
First, all ten mosaics (five for each of the 3.6/rm and 
4.5/im bands) were aligned to the same reference frame. 
We then cut them into a slightly overlapping 5x6 grid 
of 3940 X 3940 pixel subimages (with 20 pixel overlaps). 
Because the Spitzer point spread function (PSF) is “tri¬ 
angular” with diffraction spikes, we pre- flltered these im¬ 
ages to obtain Gaussian-like PSFs Isee IKozlowski et al.l 
l 2010 bl for details), which greatly helps the difference 
image analysis (DIA) software obtain clean subtrac¬ 
tions and photometry. Having the corrected images, we 
chopped each of them into even smaller 2 x 2 grids of 
1990 X 1990 pixel chunks (again with 20 pixels of over¬ 
lap) to speed up the DIA analysis. We stacked all five 
epochs in order to create template images and then use 
DIA to match and subtract the reference images from 
each individual epoch. 

We measured the fluxes of the objects in the SD¬ 
WFS catalogs on the template image to photometri- 
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TABLE 1 

The 3.6 Micron Light Curves for the DIBS Variability Catalog. 


Measurements Uncertainties 


RA 

(deg) 

Dec 

(deg) 

ei 

(mag) 

e2 

(mag) 

£3 

(mag) 

€4 

(mag) 

es 

(mag) 

ei 

(mag) 

£2 

(mag) 

£3 

(mag) 

£4 

(mag) 

es 

(mag) 

217.659222 

32.452357 

18.677 

18.742 

18.777 

18.920 

18.977 

0.141 

0.139 

0.137 

0.131 

0.153 

217.730372 

32.452616 

19.048 

18.793 

19.001 

19.217 

19.029 

0.158 

0.159 

0.155 

0.147 

0.173 

217.787424 

32.452656 

19.349 

19.573 

20.245 

19.596 

99.999 

0.207 

0.204 

0.196 

0.192 

99.999 

217.628739 

32.452706 

19.119 

19.483 

19.180 

19.391 

19.950 

0.179 

0.175 

0.174 

0.166 

0.190 

217.618530 

32.452833 

18.907 

18.853 

19.024 

18.676 

18.973 

0.147 

0.146 

0.143 

0.139 

0.160 


Note. — ei with i — 1, ...,5 are the epoch numbers. The error code for magnitudes, reflecting no measurement, is 
99.999. 

(This table is available in its entirety in a machine-readable form in the online journal. A portion is shown here for 
guidance regarding its form and content.) 


TABLE 2 

The 4.5 Micron Light Curves for the DIBS Variability Catalog. 


Measurements Uncertainties 


RA 

(deg) 

Dec 

(deg) 

(mag) 

£2 

(mag) 

£3 

(mag) 

£4 

(mag) 

es 

(mag) 

ei 

(mag) 

^2 

(mag) 

£3 

(mag) 

64 

(mag) 

^5 

(mag) 

217.659222 

32.452357 

18.875 

19.046 

19.653 

18.559 

18.317 

0.192 

0.158 

0.170 

0.172 

0.164 

217.730372 

32.452616 

18.844 

18.433 

19.255 

18.652 

18.526 

0.170 

0.143 

0.152 

0.150 

0.143 

217.787424 

32.452656 

99.999 

99.999 

99.999 

19.855 

19.689 

99.999 

99.999 

99.999 

0.321 

0.304 

217.628739 

32.452706 

18.837 

18.673 

19.173 

19.397 

19.468 

0.206 

0.171 

0.184 

0.177 

0.167 

217.618530 

32.452833 

19.314 

19.029 

18.570 

18.273 

18.591 

0.173 

0.143 

0.161 

0.158 

0.147 


Note. — ei with 2=1, ...,5 are the epoch numbers. The error code for magnitudes, reflecting no measurement, is 
99.999. 

(This table is available in its entirety in a machine-readable form in the online journal. A portion is shown here for guidance 
regarding its form and content.) 


cally align them to the SDWFS catalog. They were 
measured using the PSF fi tting routines from the DIA 
package of iWozni^ (|200Clll . DIA typically returns un¬ 
derestimated uncertainties, hence they were corrected 
with (Tnew = -\/(C X ctdia)^ + where C, = 5.1 and 
e = 0.00324 were set so the new error bars dnew match the 
dispersions for non-variable field objects. After subtrac¬ 
tion of the template, the difference fluxes at the positions 
of all sources were measured and added to the fluxes mea¬ 
sured on the template, thereby creating the light curves 
(Tables [Hand[5]). Our catalogs are limited to the objects 
detected on the template image. A transient which does 
not produce a detectable source in the (averaged) tem¬ 
plate image and is not present in the SDWFS catalog 
will not be identified. The brightest objects create mul¬ 
tiple variable artifacts in their vicinity. We masked the 
regions around these objects using the Two Mic ron All 
Sky Survey (2MASS) catalog (|Cutri et al.l 1200311 . where 
the masking region radius was 260(1 — K/11) arcsec for 
K <\1 mag stars. 

The absolute magnitudes were derived us¬ 
ing a standard ACDM cosmological model with 
{Ho,nM,^vac,^k) = (70 km s~^ Mpc~^, 0.3, 0.7, 0.0 ) 
and using K-corrections derived from lAssef et al.l (1201011 . 

3. VARIABILITY 

We were interested in sources with at least three epochs 
(Yip ^ 3) measured in both the 3.6 /rm and 4.5 /rm 
bands. From the fluxes we compute the standard de¬ 
viations u([A])i as a function of apparent magnitude 
m{[X])i (Figure (H panel (a)). We also calculated the 
variability covariance CI 2 and Pearson’s correlation r* 


coefficients between the two bands, where 

Yep 

= - {M[x]r))\ ( 1 ) 

^ j=l 

-j -Vep 

^^2 = E - (^([ 3 . 6 ])^)) 

®p ^ J=1 

X (m([4.5])} - (m([4.5])*)) , (2) 

i ^ ^12 (-O') 

u([3.6])*u([4.5])*’ ^ ’ 

and (to([A])*) is the average magnitude for the ith light 
curve in band X. A correlation coefficient of r = 1 
(—1, 0) means that the 3.6 /im and 4.5 ^m light curves 
are perfectly correlated (anti-correlated, uncorrelated). 
We expect the majority of sources in the field to be 
non-variable galaxies, hence the distribution of u([A])i 
is dominated by contributions from noise and not true 
variability. In order to select truly variable objects, we 
defined a variability significance threshold. Let Um([A]) 
be the median dispersion as a function of magnitude rrii 
and cr([A]) be the dispersion of the u([A])j values around 
the median. The more distant the source is from the me¬ 
dian, (uOAJjj — Vm{[X\)) / a{[X\) > 0, and the more its 
3.6 /iin and 4.5 /im light curves are correlated, the more 
likely it is to be the true variable. The typical rms of 
mid-IR quasar variability (marked by dots in Fig. [T|) is 
^0.1 mag, which is comparable to the uncertainties pro¬ 
vided by DIBS or SDWFS, hence we can detect only 
bright and/or significantly variable quasars. We define a 
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Fig. 1.— Panel (a): Variability in the 3.6 fim and 4.5 bands as a function of magnitude for all objects with 4 or 5 epochs of data 
shown as contours. The objects are binned into 0.1 mag and 0.05 log(u) bins. The contours are drawn for 2, 10 and 100 objects per 
bin, counting from the outer contour. The lower solid line shows the median Vm calculated in 0.1 mag bins, and the upper three dashed, 
dash-dotted and dotted lines show variability significance {y — Vm)lo' = 1, 2 and 3 (from bottom to top) with respect to this median. We 
limit our Photometry Group 2 and 3 samples to objects with [3.6] < 19.8 mag, and [4.5] < 19.3 mag (vertical dashed lines). Panels (b-d): 
Joint variability <Ji 2 as a function of [3.6] magnitude. Contours have the same meaning as in panel (a). The three horizontal lines mark the 
Variability Levels 2, 3, and 4, from bottom to top, while the vertical lines mark the mean variability dispersion at 0.05, 0.10 and 0.15 mag 
(from left to right) and the limiting magnitude (the right most solid line). In the three panels we show (b) t he AGES confirmed quasars, 
(c) X-ray detected sources, and (d) AGN candidates from fChung et ^ 1120141) . Colors in panel (d) denote the lChung et all II2014I) F-ratio 
ranges, where 10 < F < 30 is blue, 30 < F < 80 is green, and F > 80 is marked with red (see text). 


combined significance of the variances in the light curves 
as 

_2 fv{[3.6])i-vy[3.6])im^)Y , 

"“<■> = I-^(|3.6])(™,)- ) + 

^([4-5])i -t>m([4.5])(TOi) y 

a([4.5])(m.) J ^ ^ > 

which quantifies the degree to which the object deviates 
from the median variances in both bands. The best can¬ 
didates for true variables are objects with positive excess 
variances (large a 12 and v > Vm) and strong correla¬ 
tions r ~ 1 between the two bands. We define three 


Variability Levels 2, 3, and 4 (in short VL2, VL3, and 
VL4) for sources exceeding the joint variability signifi¬ 
cance CT 12 > 2, 3, and 4, respectively, that will be fre¬ 
quently used to select variable objects and discuss selec¬ 
tion yields (see FigurelU panels (b-d)). By definition half 
of the sources have negative (u([V])i — Vm{[X\))/ 
and anomalously non-variable sources could also give rise 
to high values of cri 2 . The distributions of (t!([V])i — 
nm([V])) are, however, very non-symmetric and while 
for the 3.6 /im or 4.5 /rm bands there are 7443 or 
7538 objects with (n([V])j — Vm{[X])) /(j{\X]) > 1 and 
tTi 2 (f) > 2, there are only 0 or 17 sources with (^([V])! — 
VTn{[X]))/a{[X\) < —1 and cri 2 (i) > 2. This means that 
the high variability significance sources are those with 
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correlation r 

Fig. 2. — Correlation between the 3.6 fim and 4.5 fim light curves 
for sources with measured [5.8] — [8.0] colors (“Photometry Group 
3” sources; see ai with 4 or 5 epochs of data for Variability Level 2, 
3 and 4 (cri 2 > 2, > 3 and > 4) sources. The true variables are the 
sources to the right of the vertical dashed line, in the high correla¬ 
tion region (r > 0.8). The poorly correlated (r < 0.5) sources were 
used to estimate our false positive rates inside the high correlation 
region (Tables [JJ [S] and[^. 
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Fig. 3.— Top panel: Histograms (solid lines, left scale) and cumu¬ 
lative distributions (dashed lines, right scale) for spectroscopically 
confirmed AGES quasars with data from all five epochs and vari¬ 
ability exceeding VL2 as a function of observed [3.6] magnitude. 
Red, blue, and black lines are for high correlation sources with 
r > 0.8, for r < 0.8, and all AGNs, respectively. Bottom panel: 
Ratios of histograms and cumulative distributions for highly cor¬ 
related AGNs to all AGNs. 



Fig. 4. — Mid-IR color-color distribution of the Photometry 
Group 3 sources in DIBS (green). The contours are based on 39591 
objects that are counted in 0.05 mag bins in both axes. The con¬ 
tours are at levels of 2, 10, 20, 50, 100, and 150 (counting from 
the outside). In order to understand color-color selections, we plot 
the lAssef et al.l II2010I) color-color tracks as a function of redshift 
for common object classes such as AGNs (red dotted line), spiral 
(Sbc; solid blue line), elliptical (Ell; dashed blue line), and irregu¬ 
lar (Irr; dotted blue line) galaxies. We also show a combined light 
of an AGN and its host spiral galaxy, contributing light equally 
(red solid line). The dots on each track mark redshift increment 
by 0.5 (1.0) from 2 : = 0 to 3 (to 6) for galaxies (AGNs) that start 
at 2 = 0 bulls eye. The solid black trapezoid is the modified AGN 
wedge, the dashed wedge is the old AGN wedge, and the black 
solid rectangle is the stellar box. 

positive (high) variability. 

The distribution of correlation strengths r for sources 
with (J 12 exceeding a range of variability significance 
thresholds is shown in Figure [U It is obvious that the 
true variables tend to be objects with high correlations 
r > 0.8 (dashed vertical line). We can estimate the 
contamination as follows. We assume that uncorrelated 
sources (r < 0.5) in Figure [5] represent the density of 
false positives. If the number of identified variables with 
(712 > 2 (or 3, 4) and r > 0.8 is Nryo.s^ and the num¬ 
ber of sources with ai 2 > 2 (or 3, 4) and r < 0.5 is 
-^r<o .55 then we estimate that the number of false pos¬ 
itives is F = (0.2/1.5)Vr<o.5, where 0.2/1.5 is the ratio 
of the intervals r > 0.8 and —1 < r < 0.5. The fraction 
of the identified variables that are false positives is then 
F/Nr> 0 .s- In Figure [Sj we show both the histograms 
and cumulative distributions for spectroscopically con¬ 
firmed quasars from the AGES survey with the variabil¬ 
ity exceeding VL2 as a function of the observed [3.6] 
magnitude. Bright sources with highly correlated light 
curves (r > 0.8) constitute ^100% of the AGN sample, 
but this proportion decreases to ^80% with decreasing 
source brightness. The variability properties of DIBS 
sources are presented in Table [3] 

Knowledge of the mid-IR colors is key to assign¬ 
ing objects to various classes such as AGNs, galax¬ 
ies or stars (Figure ip. We use the same definitions 
as in iKozlowski et alF(|20I0bf) . where the stellar box is 
-0.15 < [3.6] - [4.5] < 0.15 and -0.15 < [5.8] - [8.0] < 
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0.15, the old AGN wedge (iStern et al.ll2005ll is defined as 


[5.8] - 
[3.6] - 
[3.6] - 


8.0] > 0.6 mag (left side), 


4.5] > 2.5 X ([5.8’ 
4.5] > 0.2 X ([5.8' 


8.0]) — 3.5 mag (right side), 
8.0]) + 0.18 mag (bottom). 


and the modified AGN wedge is defined as 

) — 0.6 mag (left side), 

) — 3.5 mag (right side), 
) + 0.08 mag (bottom). 
The mid-IR colors depend vitally on the brightness of ob¬ 
jects, therefore we divide our sources into three classes 
that increase with our ability to characterize sources in 
the mid-IR color-color plane: 


3.6 

— 

4.5 

< 2.5 

X 

( 

5.8 

— 

8.0 

3.6 

— 

4.5 

> 2.5 

X 

( 

5.8 

— 

8.0 

3.6 

- 

4.5 

> 0.2 

X 

( 

5.8 

- 

8.0 


• Photometry Group 1 (PGl) consists of all objects 
that after masking have 4 or 5 epochs of data, with 
no constraints on either their magnitudes or mag¬ 
nitude uncertainties. There are 438270 Photome¬ 
try Group 1 objects after masking, of which 9811, 
4229 and 2472 meet the Variability Level 2, 3 and 
4 criteria. 


• Photometry Group 2 (PG2) is the subset of the 
Photometry Group 1 sources for which we can mea¬ 
sure the [3.6] —[4.5] color. These sources must sat¬ 
isfy [3.6] < 19.8 mag and [4.5] < 19.3 mag with 
uncertainties in these bands smaller than 0.1 mag. 
There are 214433 Photometry Group 2 objects, of 
which 5325, 2878 and 1851 meet the Variability 
Level 2, 3 and 4 criteria. 

• Photometry Group 3 (PG3) is the subset of Pho¬ 
tometry Group 2 objects for which we can also 
measure the [5.8] —[8.0] color. These sources are 
required to have uncertainties smaller than < 0.2 
mag for the 5.8 and 8.0 /xm bands. There are 39591 
Photometry Group 3 objects, of which 1615, 840 
and 504 meet the Variability Level 2, 3 and 4 cri¬ 
teria. 


In Table m we present the significance statistics for the 
DIBS sources as a function of Photometry Group, Vari¬ 
ability Level, and the mid-IR color-color class. The most 
variable, highly correlated (r > 0.8) sources in the field 
are AGNs and stars, peaking at ai 2 ~ 50. 

4. RESULTS 

We will now methodically study the variability of mid- 
IR sources based on their classifications. 

We will start with a general description of all vari¬ 
able objects. There are 438270 Photometry Group 1 
objects (with measurements in 4 or 5 epochs in both 
bands). There are 9811 (2.2%) sources that show vari¬ 
ability cri 2 > 2 (Variability Level 2) and 1569 (0.4%) 
objects that also show a highly significant correlation 
(r > 0.8). For the sources with well-defined [3.6] — [4.5] 
colors (214433 objects), we find 5325 Variability Level 
2 objects (2.5%) and 1206 (0.6%) with a high correla¬ 
tion. Photometry Group 3, the sources with well mea¬ 
sured [3.6] — [4.5] and [5.8] — [8.0] colors, contains 39591 
sources, of which 1615 (4.1%) are Variability Level 2 ob¬ 
jects and 870 (2.2%) are also highly correlated. Increas¬ 
ing the Photometry Group, means focusing on bright 
sources, where it is easier to measure variability, so it 
is not surprising that the variable fraction increases. For 
example, for Photometry Group 1, Variability Level 2 



Fig. 5.— Left panel: Mid-IR color-color distribution of all Pho¬ 
tometry Group 3, Variability Level 2 sources in DIBS. Red points 
mark highly correlated (r > 0.8) sources, while th e open blue 
square s mark less correlated sources (r < 0.8). The IStern et all 
1120051) AGN selectio n region is shown by th e dashed wedge, the 
modified wedge from IKozlowski et al.l H2010bl) is shown as the solid 
line trapezoid, and the stellar box is marked with a rectangle. The 
contours are based on 39591 Photometry Group 3 objects, that 
are counted in 0.05 mag bins in both axes. The contours are at 
levels of 2, 10, 20, 50, 100, and 150 (counting from the outside). 
Right panel: Histogram of the fraction of the Photometry Group 
2, Variability Level 2 objects relative to all Photometry Group 2 
objects as a function of the [3.6] — [4.5] color. 

and highly correlated (r > 0.8) sources with [3.6] < 20, 
19, 18, 17, and 16 mag, we find 1569/438270 = 0.4%, 
1500/334096 = 0.5%, 1120/175827 = 0.6%, 833/68330 = 
1.2%, and 487/19132 = 2.6% variable sources, respec¬ 
tively. 

In Figure [S] we examine the distribution of variable 
Photometry Group 3 sources in the mid-IR color-color 
plane. We only consider this group because both colors 
are needed to assign classes to objects based on their lo¬ 
cation on the color-color plane. Out of 39591 objects, 870 
(2.2%) are variable (VL2) with as little as 7% contami¬ 
nation (TableSj). Most of them (652; 75%) are AGN and 
are located in the modified AGN wedge. There are also 
143 VL2 stars (16%) and 75 (9%) objects outside these 
regions, hence likely AGNs with bright host galaxies. In¬ 
creasing the significance threshold to VL3, we are left 
with 514 highly correlated variable objects, of which 376 
(73%) are AGNs, 95 (18%) are stars, and 43 (8%) are out¬ 
side these regions. The final and the highest threshold, 
VL4, selects 304 sources, of which 209 (69%) are AGNs, 
68 (22%) are stars, and 27 (9%) are other sources. 

It is easy to see that with the increasing significance 
threshold, the fractional number of AGN decreases while 
that for stars increases. This is likely due to differences 
in the typical variability amplitudes: AGNs typical vary 
by only tenths of a magnitude, while stellar variability 
can be very spectacular, especially in the infrared (e .g., 
Miras or Long Period Varia bles: seelSm ith et a l1l2002l for 
the mid-IR variability and iSoszvhski et al.l 120091 for the 
optical variability at 0.8 ^m). We inspected the mid- 
IR color-magnitude diagram (CMD; [3.6] — [5.8] color 
vs. [3.6] magnitude) and found that (l) all the PG3 
VL2 variable objects in the stellar box (but one) are 
brighter than [3.6] < 15.5 mag and (2) the stellar lo¬ 
cus at the bright end is bi-modal, with a color difference 
of 0.15 mag, and that the variable stars occupy the red¬ 
der branch. We suspect that our typical variable star 
is an RR Lyrae variable. These show mid-IR variability 
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TABLE 3 

DIBS Variability Catalog 


Object Name RA Dec [3.6] [4.5] [5.8] [8.0] Ne 'cp.e] '^[4.5] *712 r (^12 M* 

DIBS (deg) (deg) (mag) (mag) (mag) (mag) (mag) (mag) (mag^) 


J143038.21-I-322708.5 

217.6592 

32.4523 

18.56 

18.56 

99.99 

17.03 

5 

0.13 

0.51 

-0.04 

-0.65 

0.85 

1 

J143055.29-I-322709.4 

217.7303 

32.4526 

18.89 

18.51 

17.44 

16.26 

5 

0.15 

0.33 

0.01 

0.22 

0.26 

1 

J143108.98-I-322709.6 

217.7874 

32.4526 

19.52 

20.45 

99.99 

99.99 

1 

99.99 

99.99 

99.99 

99.99 

1.39 

1 

J143030.90-I-322709.7 

217.6287 

32.4527 

19.31 

18.51 

99.99 

99.99 

5 

0.33 

0.35 

0.06 

0.48 

0.28 

1 

J143028.45-I-322710.2 

217.6185 

32.4528 

18.81 

18.50 

99.99 

99.99 

5 

0.13 

0.41 

0.02 

0.29 

0.59 

1 


Note. — Right ascension and declination (magnitudes and variability measures) are provided to four (two) decimal places in order to fit 
the table into the page, w hile in the electron ic version they are in their full form. The Vega magnitudes presented in the table are measured 
in 4 arcsec apertures (see lAshbv et £011200^ . The error code for magnitudes, reflecting no measurement, is 99.999. M* is the 2MASS mask, 
where 0 (1) means affected (unaffected) by a bright star. Ne is the number of epochs used in variability calculations. 

(This table is available in its entirety in a machine-readable form in the online journal. A portion is shown here for guidance regarding its 
form and content.) 


rest-frame wavelength (/^.m) 
3 2.5 2 1.5 


0.6 



rest-frame wavelength (/xm) 

4.3 4.0 3.5 3 2.5 2 1.5 1 


■d 


-30 - 


-25 - 


cd 


-20 - 



Fig. 6. — Spectroscopically confirmed AGNs from the AGES survey. In the left (right) panel, we show their [3.6] ([4.5]) absolute 
magnitudes as a function of redshift. The points are color coded, where the black open (filled) points mark the variability significance 
<ti 2 < 1 (cri 2 > 1), blue is VL2, green is VL3, and red is VL4. We also show lines of the constant observed magnitude and the rest-frame 
wavelength at which the observed light was emitted. 


amplitudes of ~0.3 mag (|Gavrilchenko et'aLll2014D that 
are similar to our sources, and with [3.6]abs ~ — 1 mag 
they would have distances (2-20 kpc) consistent with be¬ 
ing in the Galactic halo. The primary alternative, vari¬ 
able red giants, tend to ha ve higher amplitudes (0.5- 
1.0 mag; ISmith et al.l 1200^ and would have implausi¬ 
bly large distances (100 k pc) because they typ ically have 
[3.6]abs < —9 mag (e.g.. iMeixner et TH 2 OO 6 D . Since a 
fair fraction of variable stars are bright, we also note that 
some of their variability may not be real but rather be 
due to residuals from Gaussian filtering the images. By 
increasing the significance threshold, we lose less variable 
quasars more quickly than the more variable stars. 

In Figure m we present the dependence of variability 
of the AGES confirmed quasars on the redshift, absolute 
and observed magnitudes for both bands. The most vari¬ 
able objects seem to be those with redshifts 0.5 < z < 2 
and absolute magnitudes higher than —25. 

We have also matched our DIBS variabi lity catalog to 
the XBootes X-ray catalog of the held ([Murray et al.l 


l2005t Ikenter et aTI 120051 : [Brand et al.l l2006() . MIPS de- 
tections at 24 jum (from MA GES) and radio (FI RST; 
iBecker et al.lll99^ and WSRT; Ide Vries et al.l [200^ cat¬ 
alogs. We summarize the results in Table [5] and in Fig¬ 
ure [T] AGNs are known to have X-ray emission, so it is 
not surprising to have a signihcant number of matches 
and that they mostly lie in the AG N color region (Fig¬ 
ure [71 see also lGoriian et al]l2008[) . Out of 39591 PG3 
sources, we found 2307 matches to X-ray sources, of 
which 448 (19%) were at least as variable as Variabil¬ 
ity Level 2. By increasing the variability threshold to 
VL4, we are left with 160 (6.9%) X-ray-selected, mid-IR 
variable AGNs. Quasars are also known to have strong 
emission at mid- and far-IR. We matched our variabil¬ 
ity catalog to MIPS sources detected at 24 /im to hnd 
26280 matches. Of them, 713 (2.7%) were VL2 objects 
and 238 (0.9%) were VL4 sources. We also found nearly 
1500 matches to radio catalogs. However, only 51 (3.4%) 
objects were VL2 sources and there were only 17 (1.1%) 
VL4 objects. 
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TABLE 4 

Mid-IR Variable Objects in the Decadal IRAC Bootes Survey. 


Variability 

Level 

Variability Old AGN wedge 

<712 

Modified AGN wedge 

Stellar box 

Outside* 

All objects 

Number (fraction/false positive rate**) of variable objects: Photometry Group 1 

Total 


42689 

74751 

40343 

323176 

438270 

2 

> 2 

684 (1.6%/14%) 

803 (1.1%/21%) 

178 (0.4%/37%) 

588 (0.2%/100%) 

1569 (0.4%/57%) 

3 

> 3 

377 (0.9%/ll%) 

438 (0.6%/18%) 

113 (0.3%/21%) 

256 (0.1%/100%) 

807 (0.2%/47%) 

4 

> 4 

213 (0.5%/13%) 

242 (0.3%/20%) 

83 (0.2%/17%) 

152 (0.0%/100%) 

477 (0.1%/48%) 



Number (fraction/false positive rate**) of variable objects: Photometry Group 2 


Total 


29339 

55055 

18108 

141270 

214433 

2 

> 2 

651 (2.2%/9%) 

755 (1.4%/15%) 

164 (0.9%/19%) 

287 (0.2%/100%) 

1206 (0.6%/36%) 

3 

> 3 

371 (1.3%/8%) 

427 {0.8%/14%) 

108 (0.6%/ll%) 

167 (0.1%/100%) 

702 (0.3%/34%) 

4 

> 4 

208 (0.7%/ll%) 

235 (0.4%/17%) 

80 (0.4%/10%) 

111 (0.1%/100%) 

426 (0.2%/38%) 



Number (fraction/false positive rate**) of variable objects: Photometry Group 3 


Total 


6692 

12801 

7477 

19313 

39591 

2 

> 2 

593 (8.9%/2%) 

652 (5.1%/3%) 

143 (1.9%/11%) 

75 (0.4%/38%) 

870 (2.2%/7%) 

3 

> 3 

341 (5.1%/!%) 

376 (2.9%/3%) 

95 (1.3%/4%) 

43 (0.2%/34%) 

514 (l.3%/5%) 

4 

> 4 

192 (2.9%/2%) 

209 (1.6%/3%) 

68 (0.9%/3%) 

27 (0.1%/40%) 

304 (0.8%/6%) 


Note. — * — Column “Outside” includes objects that are both outside of the “Modified AGN wedge” and the “Stellar box”. 

** — The fraction of variable sources is a ratio of the number of variable objects above a given level of (T 12 and r > 0.8 to the total number 
of objects in a respective region. The false positive rate is estimated from the number of sources above the same level of <712 but with 
r < 0.5 (see ^and Figure!^. 



Fig. 7.— Left panel: Mid-IR color-color distribution of the Pho¬ 
tometry Group 3, Variability Level 2 sources selected as AGNs 
based on their X-ray (blue), radio (red) or 24 pm emission (greei^ 
The contours and selection regions are the same as in Figure 0 
Right panel: Histograms of the fraction of the Photometry Group 
2 X-ray, 24 pm and radio sources identified as Variability Level 2 
variables as a function of the [3.6] — [4.5] color. 

iChung et al.l (|2014D classified the B ootes source s based 
on up to 17-band SED fits using the I Asset et al.l (|201Clll 
templates and provided photometric redshifts for the 
sources lacking spectroscopy. They used the F-test to 
compare fits with and without an AGN component, 
where increasing values of F make the presence of an 
AGN more probable. Because non-AGNs greatly out¬ 
number AGNs, there will be a high false positive rate for 
samples of objects at F-test values that would strongly 
indicate the p resence of an AGN if only considering a 
single object. IGhung et al.l (I2014D sources with F > 40 
almost all lie in the AGN wedge, while those with 10 < 
F < 40 track both the wedge and some of the more ex¬ 
tended color region where X-ray sources are found. We 
have matched 3794 PG3 sources with F > 10, of which 
440 (11.6%) turned out to exceed VL2 and 137 (3.6%) to 


be VL4 sources. These highly correlated variable sources 
also have very low estimated false positive rates of ~1%. 
Of the 3794 PG3 sources, 2929 (77%) lie inside the Stern 
wedge, 3282 (87%) lie inside the modified wedge, and 
only a single object (0%) in the stellar box. 

The AGES survey provides spectroscopic redshifts and 
classifi cations for roughly 23 ,500 quasars, galaxies and 
stars (iKochanek et al.l I2012f) . We matched the DIBS 
variability catalog to the AGES spectroscopic catalogs 
and obtained 2179 matches for AGNs, 490 for stars, and 
11199 for galaxies in Photometry Group 3 (Table [6]). At 
Variability Level 2, we identified 402 (18.4%) AGNs, 18 
(3.7%) stars, and 164 (1.5%) galaxies (Figure |S]). By in¬ 
creasing the variability significance threshold to VL4, we 
were left with 145 (6.7%) AGNs, 8 (1.6%) stars, and 54 
(0.5%) galaxies. 

4.1. Comparison with SDWFS 

What are the implications of adding the fifth epoch 
(this study) to the SDWFS variability catalog? 

Consider the Stern wedge sources first. The number 
of sources in the wedge is nearly identical in both sur¬ 
veys with 6677 for SDWFS and 6692 for DIBS. There 
are 496 (7.4%) SDWFS VL2 PG3 variable objects and 
593 (8.9%) DIBS sources. DIBS increases the number of 
variable AGN candidates by nearly a hundred (20%) and 
reduces the estimated false positive rate (from 6% to 2%). 
Combining the two effects, the number of AGN with ro¬ 
bustly measured variabilty is increased from roughly 470 
to 581 or about 25%. 

Increasing the variability significance threshold to VL4 
results in 204 (3.1%) AGN in SDWFS and only 192 
(2.9%) in DIBS, and the false positive rate stays at 2%. 
On the other hand, if we keep the false positive rates from 
SDWFS (i.e., use a lower significance levels for DIBS), 
we obtain 496 (SDWFS) and 936 (DIBS) PG3 sources 
for a 6% false positive rate, or 204 (SDWFS) and 678 
(DIBS) PG3 objects for a 2% false positive rate. The 
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Fig. 8 . — Left panel: Mid-IR color-color distribution of the Pho¬ 
tometry Group 3, Variability Level 2 sources with spectroscopic 
classification from AGES. Blue, red, and green symbols correspond 
to objects with quasar, galaxy, and stellar spectroscopic class, re¬ 
spectively. The contours and selection regions are the same as in 
Figure [S] Right panel: The fraction of objects with stellar, galaxy 
or QSO template matches that are the Photometry Group 2, Vari¬ 
ability Level 2 variables as a function of the [3.6] — [4.5] color. 

number of variable AGN detections is roughly doubled 
in DIBS as compared to SDWFS at fixed false positive 
rates. 

We observe a similar situation for the modified AGN 
wedge. There are 584 (4.6%) SDWFS VL2 sources and 
652 (5.1%) for DIBS (12% more) out of 12741 and 12801 
PG3 sources, respectively. Along with the increased to¬ 
tal number of these variable AGNs in DIBS, we observe 
a decreased false positive rate from 7% to 3%. When 
increasing the variability significance threshold to VL4, 
we have 224 (1.8%) and 209 (1.6%) sources, respectively. 
For the fixed false positive rates from SDWFS we have 
584/997 (SDWFS/DIBS) sources for a 7% false positive 
rate or 224/719 for a 3% false positive rate. 

We also study the content of the stellar box, where 
there are 7567 (SDWFS) and 7477 (DIBS) objects. 
While in SDWFS we have 86 (1.1%) VL2 stars, there 
are 143 (1.9%) in DIBS. Increasing the variability sig¬ 
nificance threshold to VL4 results in a collection of 21 
(0.3%) stars in SDWFS and 68 (0.9%) in DIBS. 


5. MID-IR QUASAR STRUCTURE FUNCTIONS 

Quasa rs are well know n as aperiodically variable 
sources (|Ulrich et al.l 1199711 . Such brightness changes 
are now routinely and successfully modeled using 


201Ct IMacLeod et al.l 

20101 

Butler & Bloom 2011: 

MacLeod et al.l 120111 

Zu et al. 

2011: IMacLeod et al.l 


some departures on very short time scales and low 
amplitudes that are presently only probed by Kepler 
(|Mushotzkv et al.l 120111 iKasliwal et al.l I2015I1 . Because 
our light curves are very sparsely sampled, we are re- 
strict ed to comparing averag e structure functions (SF; 
e.g.. iKawaguchi et al'l Il998t iVanden Berk et al.l 120041 : 
Ide Vries et al.l 120051 iSchmidt et al.ll2010ll . which quan¬ 
tify the mean magnitude difference at the observed (r = 
\tj — ti\) or rest (r = \tj — + z)) frame time differ- 

ence (or “lag”) b etween two epochs. We define it as in 
iKozlowski et al.l (|2010bll . 




where Nqso and Ngai are the number of quasars and 
galaxies (providing the variance/noise estimate) used, 
and m{ti) — is the magnitude difference for a given 

time-lag r. For each QSO, we match four galaxies having 
nearly identical brightness that serve as a noise estimate. 
At the relevant magnitudes, galaxies are essentially un¬ 
resolved in the DIBS images. We calculated 1000 boot¬ 
strap resamplings for both the quasar and galaxy lists, 
and we report here the median values along with the 
1(7 (68.3%) uncertainties fTables IMTH. We do not ex¬ 
pect a strong AGN variability signal at short time lags, 
which makes these estimates very sensitive to any prob¬ 
lems with the estimated noise. Having five epochs, we 
measure n(n — l)/2 = 10 different time-lags, blurred by 
the (1 -I- z) factor in the rest frame. We model the SF by 
a power-law of the form 


5(r) = ^o(^y, (6) 

adopting tq = 2 (4) years for the rest-frame (observed- 
frame) estimates. 

The AGES survey provides us with over 23500 sources 
in the field that have measured redshifts and/or type 
classification (quasar, galaxy, or star). Because the num¬ 
ber of galaxies in this magnitude limited survey falls 
off very quickly with increasing redshift, all sources be- 
yond the redshift z > 1 are quasars (see Figure 6 in 
iKochanek et al.l 1201211 . We will frequently use this in¬ 
formation, in the selection of “true” quasars and in the 
calculation of the observed and rest-frame structure func¬ 
tions. 

For the structure function, we must analyze both the 
detected and undetected variable AGNs in DIBS. We do 
this by computing the rest-frame structure function of all 
sources in a given sample with [3.6] < 18 mag. We con¬ 
sider several cases independent of their variability level. 

5.1. SF of the AGES quasars 

We calculated the rest-frame structure functions for 
spectroscopically confirmed quasars from AGES (Ta¬ 
ble [7]). The SF slopes are steeper than in the optical, 
with 7 = 0.45, and the amplitudes are Sq = 0.12 mag 
at To = 2 years (Figure [St, i n agree ment with our earlier 
findings in iKozlowski et al.l (l2010bll . The observed frame 
structure functions show slightly shallower slopes with 
7 Rs 0.4 and amplitude S'o = 0.11 mag at tq = 4 years. 

Let us now consider the dependence of the structure 
function on the redshift (or the emission wavelength) and 
the absolute luminosity of AGES quasars. We select 
three redshift ranges of 0.7 < z < 1.2, 1.2 < z < 1.5, 
and 1.5 < z < 2.0, each containing 248-336 AGNs. The 
sample of quasars in each redshift bin was then divided 
into two equal halves after sorting them by absolute mag¬ 
nitude, hereafter called the faint and bright subsamples. 
For the bright subsample, we observe constant SF slope 
across the three redshift bins (rest frame wavelengths 
A = 1.2-2.1/im) in both bands with mean 7 « 0.44. We 
observe a mild decrease of the amplitude Sq as a function 
of increasing redshift (correlation with wavelength). 
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TABLE 5 

Mid-IR Variability of Photometrically Selected AGNs. 


Variability 

Variability 

X-ray 

FIRST 

WSRT 

Chune: et al. 120141 

Level 

<T12 




Number (fraction/false positive rate) of variable objects: Photometry Group 1 

Total 


3667 

376 

2300 

9258 

2 

> 2 

469 (12.8%/!%) 

11 (2.9%/3%) 

42 (1.8%/6%) 

454 (4.9%/3%) 

3 

> 3 

286 (7.8%/l%) 

5 (1.3%/2%) 

27 (l.2%/6%) 

254 (2.7%/2%) 

4 

> 4 

163 (4.4%/l%) 

4 (l.l%/3%) 

14 (0.6%/9%) 

139 (1.5%/2%) 



Number (fraction/false positive rate) of variable objects: Photometry Group 2 

Total 


3258 

342 

2082 

5550 

2 

> 2 

468 (14.4%/1%) 

11 (3.2%/3%) 

42 (2.0%/6%) 

451 (8.1%/2%) 

3 

> 3 

286 (8.8%/l%) 

5 (1.5%/2%) 

27 (1.3%/5%) 

254 (4.6%/l%) 

4 

> 4 

163 (5.0%/!%) 

4 (1.2%/3%) 

14 (0.7%/8%) 

139 (2.5%/l%) 



Number (fraction/false positive rate) of variable objects: Photometry Group 3 

Total 


2307 

206 

1278 

3794 

2 

> 2 

448 (19.4%/1%) 

11 (5.3%/2%) 

40 (3.1%/!%) 

440 (11.6%/1%) 

3 

> 3 

277 (12.0%/0%) 

5 (2.4%/2%) 

25 (2.0%/2%) 

250 (6.6%/l%) 

4 

> 4 

160 (6.9%/0%) 

4 (1.9%/3%) 

13 (1.0%/4%) 

137 (3.6%/l%) 


Note. — The fraction of variable sources is a ratio of the number of variable objects above a given level of (T 12 
and r > 0.8 to the total number of objects in a respective class. The false positive rate is estimated from the 
number of sources above the same level of < 7^2 but with r < 0.5 (see ^and Figure!^. 


TABLE 6 

Mid-IR Variability by AGES Spectroscopic Classification. 


Variability 

Variability 

QSOs 

Stars 

Galaxies 

Level 

cri2 




Number (fraction/false positive rate) variable objects: Photometry Group 1 

Total 


2911 

1102 

18039 

2 

> 2 

415 (14.3%/!%) 

22 (2.0%/16%) 

181 (1.0%/16%) 

3 

> 3 

245 (8.4%/l%) 

13 (1.2%/16%) 

110 (0.6%/15%) 

4 

> 4 

146 (5.0%/!%) 

9 (0.8%/10%) 

62 (0.3%/21%) 



Number (fraction/false positive rate) of variable objects: Photometry Group 2 

Total 


2793 

960 

17368 

2 

> 2 

415 (14.9%/!%) 

22 (2.3%/ll%) 

178 (1.0%/13%) 

3 

> 3 

245 (8.8%/l%) 

13 (1.4%/10%) 

109 (0.6%/14%) 

4 

> 4 

146 (5.2%/l%) 

9 (0.9%/8%) 

62 (0.4%/20%) 



Number (fraction/false positive rate) of variable objects: Photometry Group 3 

Total 


2179 

490 

11199 

2 

> 2 

402 (18.4%/1%) 

18 (3.7%/5%) 

164 (1.5%/7%) 

3 

> 3 

240 (ll.0%/l%) 

11 (2.2%/4%) 

100 (0.9%/6%) 

4 

> 4 

145 (6.7%/l%) 

8 (1.6%/3%) 

54 (0.5%/9%) 


Note. — The fraction of variable sources is a ratio of the number of variable objects above a given level of (T 12 and r > 0.8 
to the total number of objects in a respective class. The false positive rate is estimated from the number of sources above the 
same level of <712 but with r < 0.5 (see m and Figure O. 


For the faint subsample, the SF slope is most likely con¬ 
stant at 7 ~ 0.46, but we observe a decline to 7 ~ 0.37 in 
both bands in the high-z bin. From Figure IHl we see that 
with increasing redshift we are losing intrinsically faint 
AGNs (that seem to be more variable) due to the ob¬ 
servational magnitude limit. The observed decline of SF 
slope in the faint sample may simply reflect this fact, be¬ 
cause we do not observe such a decline in the more com¬ 
plete bright subsamples. For both the bright and faint 
subsamples, the amplitudes Sq decline to higher redshifts 
and decreasing emission wavelengths. Also fainter sub¬ 
samples show larger amplitudes than for the bright sam¬ 


ples. 

We repeated the SF calculation in equal magnitude- 
redshift bins by dividing the AGES quasars into abso¬ 
lute magnitude bins (spaced every 1 mag) and redshift 
bins (spaced every logz = 0 . 2 ), and considered only 
bins with a minimum of 10 AGNs. The amplitude at 
a fixed time scale increases with increasing rest frame 
wavelength (or decreasing redshift). The slope is either 
constant with the rest frame wavelength (or redshift) or 
it is very weakly correlated with the rest frame wave¬ 
length (weakly anti-correlated with redshift). 

To increase the statistical sample, we also define an- 
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rest-frame time-lag (years) 


Fig. 9.— Rest-frame mid-IR AGN structure functions for AGES quasars that are brighter than [3.6] < 18 mag. The solid lines repr esent 
the fit to the points. For comparison we also show the mid-IR stru cture functions for SDWF S (dotted line: IKozlowski et aLll^lObl) and 
the optical f-band structure function for SDSS quasars (dashed line: IVanden Berk et al.1l2004l) . 



Fig. 10.— Rest-frame mid-IR AGN structure functions for the 0.5 < z < 2 AGES quasars brighter than [3.6] < 18 mag. We divide the 
sample of 1000 AGN in half after sorting them by their absolute magnitudes and calculate SF for the brighter half (red) and the fainter 
half (blue). It is clear that fainter AGN show higher amount of variability then the brighter ones. 


other larger sample of 1000 AGNs with 0.5 < z < 2. We 
again sort it by absolute brightness and split it into two 
brightness bins (Figure [10]). We do not find the slopes of 
the fainter AGNs to differ from the brighter ones (Fig¬ 
ure [13, but the amplitudes are larger for fainter AGNs. 
This sample consists of 40% highly correlated objects 
and 60% with low or no correlation. For the highly cor¬ 
related objects, we again observe larger amplitudes for 
fainter AGNs, but also steeper SF slopes than for the 
bright AGNs. 

5.2. SF of Mid-IR-selected Quasars 

We can also study the structure function of all AGN 
candidates lying inside the modified AGN color wedge. 
These sources must meet the PG3 requirement but again 
with no restriction on their variability. There are 999 
sources inside the wedge that are brighter than [3.6] < 


18 mag and at redshifts z > 1. The latter requirement 
guarantees that we are dealing with true AGNs, because 
galaxies are too faint to be observed at such redshifts in 
AGES. The rest frame structure function slopes are 7 = 
0.43-0.48 and amplitudes of So = 0.12 mag (TablejS]), in 
agreement with the true AGES quasars. 

We are also interested in the observed structure func¬ 
tions for [3.6] < 18 mag PG3 sources inside the AGN 
color wedge with either z > 1 (999 sources) or ignor¬ 
ing redshift information (11302 sources). The observed 
slopes for the first sample are 7 = 0.42 ± 0.08 and 
7 = 0.39 ± 0.04 for the 3.6 /im and 4.5 /im bands. 
The amplitudes Sq at the time scale to = 4 years are 
So = 0.102 ± 0.012 mag and Sq = 0.107 ± 0.008 mag, 
respectively. The second sample is characterized by 
SFs with 7 = 0.37 ± 0.12 and 7 = 0.34 ± 0.08 for 
the 3.6 ^m and 4.5 /rm bands, and the amplitudes of 
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So = 0.064 ± 0.012 mag and S'o = 0.064 ± 0.009 mag, 
respectively. 

5.3. SF of X-ray-selected Quasars 

We analyzed the X-ray-matched sources with known 
redshifts from AGES that lie inside the modified AGN 
wedge. We calculated the rest-frame structure functions 
(TableO and obtained the following slopes 7 ~ 0.46 and 
amplitudes S'o = 0.13 mag for AGNs at z > 1, and ob¬ 
tained slopes 7 ~ 0.44 and amplitudes So = 0.14 mag for 
all AGNs. These values again stay in general agreement 
with the measurements for the confirmed AGES quasars. 

The observed structure functions for all [3.6] < 18 mag 
PG3 X-ray sources inside the AGN color wedge, with 
or without known redshifts (1720 sources) have slopes 
7 = 0.39 ± 0.07 and 7 = 0.37 ± 0.04 for the 3.6 ^m 
and 4.5 /rm bands, and the amplitudes Sq at the time 
scale to = 4 years are So = 0.109 ± 0.011 mag and Sq = 
0.122 ± 0.008 mag, respectively. 

5.4. SF of Radio-seleeted Quasars 

We also inspected the variability of 183 radio sources 
from inside the modified AGN wedge with known red¬ 
shifts. We measure the rest-frame slopes 7 = 0.39±0.10 
for 3.6 fim and 7 = 0.36±0.08 for 4.5 fim bands. This se¬ 
lection returns slightly shallower slopes than for the con¬ 
firmed AGNs. The amplitudes at 2 years are So = 0.08 
and So = 0.09 mag, respectively. The observed struc¬ 
ture functions for 555 sources from inside the modified 
AGN wedge with either known or unknown redshifts have 
slopes 7 = 0.32±0.08 and 7 = 0.36±0.09 and amplitudes 
at the time scale to = 4 years So = 0.060±0.008 mag and 
So = 0.058 ± 0.008 mag, in the two bands, respectively. 

5.5. SF of Quasar-candidates from \Chuna et al\ 11201 A ) 

iGhung et al.l (j2014t l provide a catalog of AGN can- 
didates bas e d on up to 17-band SED fits using the 
lAssef et ’HI (|2010t l templates. We have kept only the 
sources with [3.6] < 18 mag, F-ratio F > 10 or F > 30 
for which we adopted the AGES spectroscopic redshifts 
and if unavailable then we used the estimated photomet¬ 
ric redshifts. We find the rest-frame slopes 7 ss 0.37 and 
amplitudes So ~ 0.08 mag for the F > 10 sources, and 
7 « 0.40 and amplitudes So ~ 0.09 mag for the F > 30 
sources. These values are slightly lower than the ones 
for the confirmed AGNs. This is most li kely an indi¬ 
cation of the level of contamination in the iChung et al.l 
(|2014li selection method. In fact, we can use the am- 
plitud e of the structure f unction to estimate the purity 
of the IGhung et al.l (|2014ll quasar candidate sample as a 
function of the F-test values. We assume that the spec¬ 
troscopically confirmed AGES AGN are representative 
of AGN variability with Strue = 0.12 mag (0.13) mag at 
to = 2 years and 7 = 0.47 (0.45) in the 3.6 ^m (4.5^m) 
band. 

If we measure amplitude S{F) (at fixed slope) for a 
sample where fraction / is composed of sources varying 
with an average amplitude of Ftrue and fraction 1 — / are 
non-variable, then the fraction of real sources is /(F) = 
S'(F)^/S'/^^g. The scaling of /(F) with F is robust, but 
the quantitative estimate depends strongly on using the 
correct value of Ftrue, particularly as S'(F) —>■ Ftrue- 
The formal uncertainty is u/ = 2/(Ttrue/<S'true so even 


TABLE 7 

Mid-IR Structure Functions for confirmed 
AGES Quasars ([3.6] < 18 mag) 


Band 

Structure Function Parameters 
7 S'o (mag) 

Nobj 

rest-frame (tq = 2 yrs) 

3.6 fim 

0.468 ± 0.073 

0.118 ±0.007 

1518 

4.5 fim 

0.445 ± 0.042 

0.128 ±0.004 

1518 

3.6 fim 

fixed 0.5 

0.113 ±0.006 

1518 

4.5 fim 

fixed 0.5 

0.134 ±0.006 

1518 

observed-frame (tq = 4 yrs) 

3.6 fim 

0.373 ± 0.081 

0.110 ±0.014 

1518 

4.5 fim 

0.423 ±0.041 

0.121 ±0.007 

1518 


our roughly (Ttrue/*S'true — 5% measurement uncertainties 
represent a 10% un certainty in / ( F) as /(F) —>• 1. 

If we consider the IGhung et al.l (|2014D sources in “log¬ 
arithmic” bins of F = 1 — 3,3 — 10, 10 — 30 and > 30, we 
find S'(F) = 0.025, 0.029, 0.070 and 0.091, respectively, 
at 3.6/i,m. This then implies / = 0.04 ±0.01, 0.06 ±0.01, 
0.35 ± 0.05 and 0.59 ± 0.04. The results for 4.5 /im are 
similar. We clearly see that the purity of the sample 
increases with F. The samples with F < 10 are domi¬ 
nated by contaminants, as expected. While these values 
of F would imply a significant AGN component when 
considering a single object, they are not as significant 
when used for a sample of objects where only ^ 1 % of 
the s ources are AGN (see the discussion in IGhung et al.l 
I2014D . The low fraction found for F > 30 sources sug¬ 
gests that the uncertainties in Ftrue are dominated by 
systematic effects (e.g., are the AGES AGN representa¬ 
tive?) rather than the measurement errors. The fraction 
/(F) is, of course, a crude statistical estimate and likely 
subject to a number of systematic effects. 

5.6. Optical vs. Mid-IR SFs 

From studies of optical quasar variability, we know that 
the structure function sl opes are shallower at 7 ~ 0.3 
l|Vanden Berk et ^120041) than the mid-IR ones reported 
here. In Figure [U we compare the optical and mid-IR 
structure functions. It is clear that at short time scales 
the variability amplitude in the mid-IR is significantly 
lower than that in optical, and the two match each other 
only at time-lags of several years. This is consistent with 
the wavelength depen dent smoothing seen in NGC 2617 
(jShappee et al.ll2014D . 

6. SUMMARY 

In this paper, we have combined the images from SD- 
WFS with a fifth epoch from DIBS, extending the time 
baseline of the survey to 10 years. We performed DIA 
photometry on the five epochs of DIBS data, providing 
the light curves and variability measures for half a million 
mostly extragalactic sources. 

We cross-correlated the DIBS variability data with the 
AGES redshift survey. X-ray, MIPS, radio data, and 
photometrically selected AGN candidates from the same 
area of the sky. In comparison to SDWFS, the variabil¬ 
ity yields are by 15-25% larger for X-ray and 24 p,m 
quasars, and nearly identical for radio sources exceeding 
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TABLE 8 

Mid-IR Structure Functions for AGN 

CANDIDATES INSIDE THE MODIFIED AGN WEDGE 
(sources WITH Z > 1 AND [3.6] < 18 MAG) 


Band 

Structure Function Parameters 

7 So (mag) 

Nobj 

rest-frame (ro = 2 yrs) 

3.6 fim 

0.479 ± 0.078 

0.116 ±0.007 

999 

4.5 fim 

0.425 ± 0.043 

0.122 ±0.004 

999 

3.6 fim 

fixed 0.5 

0.111 ±0.005 

999 

4.5 fim 

fixed 0.5 

0.122 ±0.003 

999 

observed-frame (tq = 4 yrs) 

3.6 fim 

0.423 ± 0.080 

0.102 ±0.012 

999 

4.5 fim 

0.392 ±0.042 

0.107 ±0.008 

999 


TABLE 9 

Mid-IR Structure Functions for X-ray-detected 
sources inside the modified AGN WEDGE 
([3.6] < 18 mag) 


Structure Function Parameters N 


Band 

z > 1 

7 

So (mag) 




rest-frame (tq = 

2 yrs) 


3.6 fim 

Y 

0.463 ± 0.093 

0.125 ±0.009 

625 

4.5 ^m 

Y 

0.458 ± 0.042 

0.126 ±0.005 

625 

3.6 ^m 

Y 

fixed 0.5 

0.128 ±0.008 

625 

4.5 fim 

Y 

fixed 0.5 

0.130 ±0.010 

625 

3.6 fim 

N 

0.435 ± 0.082 

0.125 ±0.008 

1079 

4.5 fim 

N 

0.432 ± 0.045 

0.136 ±0.005 

1079 

3.6 fim 

N 

fixed 0.5 

0.125 ±0.007 

1079 

4.5 fim 

N 

fixed 0.5 

0.138 ±0.006 

1079 


observed-frame (tq 

= 4 yrs) 


3.6 fim 

Y 

0.412 ±0.058 

0.112 ±0.012 

625 

4.5 fim 

Y 

0.430 ± 0.047 

0.111 ±0.009 

625 

3.6 fim 

N 

0.385 ± 0.068 

0.109 ±0.011 

1720 

4.5 fim 

N 

0.372 ± 0.039 

0.122 ±0.008 

1720 


the variability significance cri 2 > 2. The yield of variable 
AGNs selected from the AGES survey seem to be also 
increased by '^20% as compared to SDWFS, while the 
number of variable stars and galaxies stay the same. The 
mid-IR variability selection detects 18% of the known 
AGNs from the AGES survey. Keeping the false pos¬ 
itive rates as in SDWFS (i.e., lowering the variability 
significance level used for the DIBS sample) doubles the 
number of selected variable AGNs. The addition of the 
fifth epoch significantly lowered our estimated false pos¬ 
itive rates. The primary reason for the low overall yields 
is that quasar variability amplitudes are comparable (for 
[3.6] < 18 mag) to the photometric quality of the DIBS 
data. 

The main motivation of this paper was to study the 
variability parameters (the SF slope and amplitude) as 
a function of physical quantities (luminosity or wave¬ 
length) or selection method (true AGNs, or mid-IR, X- 
ray, or radio selected, or from SED fitting). The variabil¬ 
ity analysis of ~1500 spectroscopically confirmed sources 
from the AGES survey returned the rest-frame structure 
function slopes in both bands of 7 ~ 0.45 and ampli- 


TABLE 10 


Mid-IR Structure Functions for 
[3.6] < 18 mag AGN candida tes 



fromIChung 

ET AL.I (120141') 



Structure Function Parameters 

N 

Band 

7 

So (mag) 


F > 30, rest-frame (tq = 2 yrs) 

3.6 fim 

0.431 ±0.081 

0.087 ± 0.007 

1518 

4.5 fim 

0.384 ± 0.077 

0.106 ±0.008 

1518 

3.6 ^m 

fixed 0.5 

0.084 ± 0.005 

1518 

4.5 ^m 

fixed 0.5 

0.100 ±0.005 

1518 


F > 30, observed-frame (tq = 4 yrs) 


3.6 fim 

0.395 ±0.050 

0.086 ± 0.008 

1518 

4.5 ^m 

0.381 ± 0.048 

0.096 ± 0.008 

1518 


F > 10, rest-frame (tq = 2 yrs) 


3.6 fim 

0.370 ±0.072 

0.067 ± 0.006 

4160 

4.5 ^m 

0.335 ± 0.044 

0.090 ± 0.003 

4160 

3.6 ^m 

fixed 0.5 

0.075 ± 0.005 

4160 

4.5 fim 

fixed 0.5 

0.095 ± 0.006 

4160 

F > 10, observed-frame (tq = 4 yrs) 

3.6 fim 

0.364 ± 0.059 

0.069 ± 0.007 

4160 

4.5 fim 

0.382 ±0.056 

0.078 ± 0.008 

4160 


tudes So = 0.11 mag (at tq = 2 years). Except for the 
small sample of radio-selected AGNs the mid-IR struc¬ 
ture functions of AGNs are significantly steeper than in 
optical (7 ~ 0.3). The optical and mid-IR variability am¬ 
plitudes match each other for time-lags of several years, 
but for the shorter lags the mid-IR variability amplitudes 
are significantly lower than the ones in optical. The lack 
of short time scale variability can be interpreted as the 
accretion disk or dust torus smoothing out short-time 
scale variations. This is easily understood, as the accre¬ 
tion disk radius is ~ 4 x (Arest/ttm)"'/^ light days for 
a typical black hole mass of 10 ® solar masses, radiative 
efficiency of 0.1, and Eddington ratio of 1/3, hence the 
optical variations originate in a region of 1 light-day in 
size and the mid-IR ones from a region 25 times bigger. 

We also explored the dependence of the structure func¬ 
tion slope and amplitude as a function of the luminos¬ 
ity and rest-frame wavelength. We find significant anti¬ 
correlation (lack of correlation) of the AGN amplitude 
(slope) with luminosity. In accordance with optical stud¬ 
ies, the brighter the AGN the lower its variability ampli¬ 
tude. We also find no clear evidence for the correlation of 
the variability slope with the rest-frame wavelength, but 
the range of rest-frame wavelengths explored is rather 
narrow (A = I.2-2.I/xm). On the other hand, there is 
a correlation of the SF amplitude with the rest-frame 
wavelength. 
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